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Abstract. DNA Self-Assembly has emerged as an interdisciplinary field 
with many intriguing applications such DNA bio-sensor, DNA circuits, 
DNA storage, drug delivery etc. Tile assembly model of DNA has been 
studied for various computational primitives such as addition, subtrac¬ 
tion, multiplication, and division. Xuncai et. al. gave computational DNA 
tiles to perform division of a number but the output had integer quo¬ 
tient In this work, we simply modify their method of division to im¬ 
prove its compatibility with further computation and this modification 
has found its application in computing rational numbers, both recurring 
and terminating, with computational tile complexity of 0(l) and 0(/i) 
respectively. Additionally, we also propose a method to compute square- 
root of a number with computational tile complexity of 0(n) for an n 
bit number. Finally, after combining tiles of division and square-root, 
we propose a simple way to compute the ubiquitously used irrational 
number, tt, using its GregoryLeibniz infinite series. 


1 INTRODUCTION 

In the last two decades, DNA has been used for archival data storage [2], molecu¬ 
lar logic circuits [3], building complex 2D structures from single stranded DNA [4] 
as well as synthetic 3D polyhedra [5]. This high utility of DNA is because of 
its known nano-scale structure, programmable nature, and flexibility [6]. Using 
double crossover (DX) and triple crossover (TX) molecules, Winfree showed that 
DNA self-assembly is Turing Universal la mini Ho]. To approximate self-assembly 
of DNA, he also proposed two models - abstract Tile Assembly Model (aTAM) 
and kinetic Tile Assembly Model (kTAM). Many well-known problems have been 
solved using tile assembly model muniiniiTiiTS]. In particular, tile-assembly 
models to perform arithmetic computations such as addition, multiplication and 
factor were shown by Yuri Burn [13] [TU [15] . Addition and multiplications op¬ 
erations had computational tile complexity of 0(l) and running complexity of 
0(n) [g where as the factoring numbers nondeterministically had 0(l) dis¬ 
tinct components m Similarly, Xuncai et. al. also came up with computational 
tiles for subtraction and division of integers, again, using 0(l) computational 


tiles. The running time complexity is 0(n) for subtraction and O(n^) for divi¬ 
sion [1]. Additionally, open source softwares such as XTile (it can convert any 
computational tile formula to a .tile file supported XGrow [E]), ISU TAS m 
and XtileMod (it can generate required computational tiles for performing arith¬ 
metic operations as well as the corresponding .tiles file for DNA tiles) have been 
developed mnni. 

In the present paper, we propose a method to compute rational number, ir¬ 
rational number( tt), and square-root of a number. Additionally, we also propose 
generalized method to compute division by extending Xuncai’s method. In [T], 
the output of the division was in the form of quotient and remainder which, how¬ 
ever, we have changed to decimal form. This paper could be useful for creating 
DNA based calculator. 

The paper is organized as follows. In Section 2, we briefly introduce aTAM. 
Section 3 provides basic computational tiles of Xuncai et. al. division. This tile 
systems are also used in computing square-root, rational number and tt with 
some minor modifications, in subsequent sections. In Section 4, we provide com¬ 
putational tiles to perform division in rational form. In Section 5, we have given 
an interesting additional tile system, the insertion tile system, required to com¬ 
pute square-root. Also, we show how to compute square-root. In Section 6, we 
have given method to compute rational number and tt as application of the 
square-root tile system. In Section 7, we briefly tell about our command line 
tool which can be used to generate .tile file for computing square-root of a num¬ 
ber. Section 8 concludes paper with general remark. 

2 BACKGROUND 

In [7], Winfree gave abstract Tile Assembly Model (aTAM) and kinetic Tile 
Assembly Model (kTAM) to understand the relation between self-assembly and 
computation. The basic building block of aTAM is a unit square tile. The edges 
of these tiles are nothing but sticky ends of DNA which can combine with other 
sticky ends of same edge label or edge color. Every tile has 4 edges - North, 
South, East and West - with a corresponding edge color and glue strength. Also, 
it is to be noted that, a tile cannot be rotated. In this paper, most of the symbols 
used have identical meaning to those used in [T] . Mathematically, a tile system is 
represented as Sr = (T^S^g^r)^ as shown in Eig. where T is the tile system, 
S is the seed configuration of this tile assembly, g is the glue strength of an 
edge, and r (always > 0) is the threshold temperature. The tile assembly starts 
growing from seed tile and input tiles at the bottom of the tile assembly system. 
A new tile will attach to the current tile if : 

— Edges have same glue color or glue label. 

— The combined glue strength of both the tiles is greater than or equal to the 

threshold temperature (r). 

Eor a given a seed configuration, we can represent above mentioned conditions 
in mathematical form. A tile t can attach at point (x, y) if following condition 
is satisfied [T]: 






S= Seed Configuration 
T= Tile System = {T1,T2,T3} 
E =Edge Colour = {x,y,p,q,r,r 



i.n.u.v} 


Fig. 1: The image shows basic concept for self-assembly of tiles. 


'^gibdDit),bdD-iiSR{x,y))) > r (1) 

deD 

Here, D indicates set of four directions i.e {N, E, S, W} and D~^ indicates 
direction opposite to D. bcinit) indicates binding domain of tile t in direction D. 
g(a, b) indicates combined glue strength of glue colors a and b. Given this set 
of condition, a tile assembly will start growing from initial tile set configuration. 
Boundary tiles with one edge of glue strength zero are used to avoid any further 

n —1 

growth of tiles. We represent a number a = by aoaia 2 a 3 ... a^_i where 

i=0 

ao is the most significant bit and a^-i is the least significant bit. For more 
details, the reader should refer to mm- 

3 TILE SYSTEMS 

Square-root is function which is provided, almost, by every existing calculator. 
However, in order to compute square-root using DNA tiles, we propose to use 
four different kinds of tiles systems, which are similar to those provided by [T], 
however, with some modifications. The similar tile systems will also be used to 
compute rational numbers and tt. 

1. Comparator Tile System 

2. Shift Tile System 

3. Shift and Subtract Tile System 

3.1 COMPARE TILE SYSTEM 

Compare tile system [1] helps us to identify relationship between two numbers 
a and b that is whether a is {>, < or =} than b. As shown in the Fig. a basic 
comparison tile consists of 12 different tiles where R = {<, > or =} depending 
on the relationship between aoaia 2 ... ak-i and 60^1 • • • For example, if a = 
IOIII 2 and b = IIOOO 2 then, as shown in Fig.|^ we take input tiles, SO, and CL 



















in seed configuration which grows to unique configuration as shown in Fig. 
Here, since uq = bo so R = {=} for first bit. For second bit, ai < bi so R = 
{<}. Once we get a {> or <}, we continue this sign till the end of system since 
irrespective of coming bits since the first number is already proved {> or <} 
than other. This proved that a < b so all the succeeding tiles will have same R. 
The proof that compare tile system is a step configuration, total k +1 steps so 
0(/c), and producer of unique tile system for the given seed configuration has 
already been given in [1]. 



Fig. 2: Comparator System: Left Tile is the general tile representation for all the 
12 comparator tiles. Here a, b, c e{0,1}. 



Fig. 3: The image shows how compare tile system grows from seed configuration 
for a = IOIII 2 and b = IIOOO 2 . 

















































3.2 SHIFT TILE SYSTEM 


Shift tile system [T] helps us to right shift the bits of a k bit number by one 
place while also padding a 0 in most significant bit. As shown in Fig. the 
system has input bit a, b, and c coming from the south end of tile. Here, a and 
b represent the input bits and c represents the right bit of previous tile (left 
tile). For example, as shown in Fig. for input a = IOIII 2 , seed configuration 
includes SO, SL, and input tiles. This seed configuration gives unique final output 
a = 0101 12 , as shown in Fig. where a is shifted right by one bit. The proof 
for the shift tile system is given in [T] as the Duplicator Tile System. This seed 
configuration also produces unique final output and is a step configuration with 
time complexity 0(/c) for a k bit input number. 
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Fig. 4: Shift System: Left Tile is the general tile representation for all the 8 shift 
tiles. 
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Fig. 5: The image shows how shift tile system grows from seed configuration for 
a=1011l2. The number a, which is to be shifted, is right bit of input tiles. 





























































































3.3 SUBTRACT SHIFT TILE SYSTEM 

Subtract & shift tile system [1] is used to subtract right input bits (divisor) 
from left input bits (divident) and then shift the right bits by one. Fig. shows 
computational tiles for subtracting and shifting a number. Here, a, b, c, and d 
are inputs where b is current right bit and c is the right bit of previous tile. 
This c bit, therefore, helps to right shift the number and xor operation helps in 
subtraction. For example, as shown in Fig. let a=10111 and b=01100. Hence, 
growth of this tile system gives a - b = 01011 as left bits and right sift version of 
b as right bits. Again, the proof and details of this method are provided in [T], 
so we will not go in to details here. But this tile system is also step configuration 
with time complexity of 0[k) taking total of /c + 1 steps for a k bit input number. 



Fig. 6: Shift and Subtract System: Left Tile is the general tile representation for 
all the 16 shift tiles. 
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Fig. 7: The growth of subtract and shift tile system from seed configuration for 
a=10111 and b=01100. 



























































































4 EXTENDING DIVISION 


Earlier, Xuncai et. al. had proposed a method to compute division using DNA 
tiles [1]. However, their idea of division considers dividing dividend by divisor to 
get quotient and remainder. Instead, here we propose a method to get output 
in different form. We give quotient in decimal form so that one can use this 
decimal quotient, subsequently, in doing further calculations using tile assembly. 
Computing division is divided in two cases as shown in the Fig. 



Fig. 8: Flow chart showing both the cases for computing division of a number. 


We have to see whether divisor is greater than dividend or vice versa. If the 
divisor is greater than dividend then we write dividend bits from its right end 
and prepend zero’s, as shown in Fig. In the other case, we write divisor bits 
from its left end and append zero’s for the remaining bits, as shown in Fig.|^ In 
the paper, we give an example of division where quotient turns out to be infinite. 
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Fig. 9: Divident = 6 and divisor = 1. Fig. 10: Divident = 1 and divisor = 6. 


Fig. p4] shows the tile assembly growth for division of 23.5 by 6. The answer 
to this division is 3.91666666... - a recurring rational number. Suppose we want 
first eight seven bits of the number. In that case, we have to give eight bit input 
to the tile assembly so that we get 5 bits before decimal place and 3 bits after 
decimal place. Note that this division is an example of case 1 as depicted in Fig.[^ 
but the number of bits are, inadvertently, same for both divisor and dividend. 





















Fig. 11: Division of a = 23.5 by b = 6 is calculated by this self-assembly. Here a is 
encoded in binary as IOIII.I 2 and b is encoded as IIO 2 . We get II.IIIO 2 as the 
output at the left side of the tile assembly model. Yellow tiles simply indicates 
decimal point. Grey tiles indicate values after decimal point. 


Note that the dot tiles used in the division process are optional. They are 
just for better understanding of the reader. Also, the square dots in the Fig. 
indicate that the assembly can be extended further if we want more answer bits. 
















































































































































The SSL and SL bits represent our answer - quotient - when seen in bottom- 
up fashion. In this case, we connect 0111110 .. .2 from the left boundary tiles and 
we place a decimal point after three bits which gives us OII.IIIO 2 = 3.916. After 
we get eight bits of quotient, we can terminate growth of this tile assembly with 
the help of boundary tiles. To sum up, we have not used any additional tiles to 
compute quotient in decimal form that is these are the exact same tile used in 
[T]. Thus, the computation uses 0(l) computational tiles and this computation 
takes 0[k‘^) steps for a k bit number. However, if we wish to terminate the 
computation then we will require 0(h) boundary tiles as shown in Fig. |ll] where 
h is the height of assembly. 


5 SQUARE-ROOT 

In this section, by combining the tile systems mentioned in Section 4, we il¬ 
lustrate how to find square-root a number. We will describe a new tile system 
required to compute square-root of a number. 


5.1 INSERT BIT TILE SYSTEM 


This tile system is used to insert a bit in a number. Fig. shows the computa¬ 
tional tiles required to do this operation. As shown in the Fig.[^ a, 6 , c, and #i 
are the inputs. The values on the east side are output values so this tile system 
grows from west to east. Here, we subtract the value of i unless it becomes zero 
and as soon at it does, we stop. We do this since we want to insert bit at 
position in the number. So, as soon as i becomes zero, we place the value of 
c there and after this tile, we shift all the remaining bits by one to right. For 
example, as shown in Fig. 13, to insert 1 at 5^^ position in IOOIOO 2 , we take 7^41 


as the east input of IL which is reduced unless 4 becomes 0 to insert 1 as right 
bit at that position and shift later bits by one place. 

Let us take an example number n = IOOIOO 2 . Suppose, we want to insert 
number 1 at fifth position. Therefore, as shown in Fig. 13, bdN{S{—k,—l)) = 


( 7 ^ 4 ,1). Now, in the second level, we keep on reducing the value unless we 
get zero. As soon as we get a zero, we insert the c bit at that position. After 
this point, we simply right shift the remaining values. 


Theorem 1 Let be an alphabet defined as 

E = {00,01,10,11,/} u 

{#ic|l < i < n — 1 and i = 2k + 1 ; e e {0, 1}} 

and Ti be the set of tiles defined over Y- let g = 1, r = 2, and S 

be a seed eonfiguration as defined in Fig. Let q be the number to insert the 
bit and let m and n be the sizes, in bits, of p and q, respeetively. If m < n, the 
number q needs to be padded to be m bits long with extra 0 in the qs high bit. Let 
Sr = (Ti,5',^,r). Then, there exists some {xo,yo)eZ 2 sueh that: S{xo + l,^o “ 
1) = Eo,S{xo - k,yo - 1) = So,S{xo - k,yo) = IL,S{xq + l,^o) = IR; for all 
0 < i < k — l,bd]sf{S{xo — i,yo — 1 )) = piqi. For all other positions (x, y), S(x, 





Fig. 12: Insert Tile System: General Tiles for inserting a bit at position from 
left side. Here a, b, c e{0,1} and 1 < i < n — 1. 



Fig. 13: The assembly inserts 1 at the 5^^ position from left in a = IOOIOO 2 . 


y) = empty. Then, S produees a unique final eonfiguration F on S to insert a bit 
at position in q. The assembly eonfiguration time is O (m). 

Proof. Consider the tile system Sr. Let input tiles be Ti = {Eq = {e, null, null, null}, 
Doo = {00, null, null, null}, IL = {null, fi^ix,$ix, null}, D 01 = {01, null, null, null}, 

Dio = {10, null, null, null}. Dll = {11^ null, null, null}. So = {$ix, null, null, null}}. 

Let the seed configuration S \ Z 2 ^ Ti be: 

Do = Sr{1, -1) 

So = SR{-k,-l) 

< IL = SR{-k,0) (2) 

m{0,l,..,k-l},SR{-i,-l)=Dp, 

For all other (x, y)eZ 2 , SR{x,y) = empty 

For the given seed configuration S as shown in Fig.[^ there is only one position 
where the incoming tile can attach. Also, given that r = 2, by induction we can 
say that Vte Ti, (bds{t), bdw{t)) is unique. Therefore, it follows that Sr produces 
unique final configuration F and it takes k -\-l steps to assemble. 

As shown in the Figj^ value of #i indicates whether the bit needs to be 
inserted or simply copied to the other side. If (#i, c) has i > 0 then the output 


























































on east will be (#(i - 1), c). If, on the other hand, we have (#0, c) as the input 
on west side, then we will replace (a, b) by (a, c) and the output on east side 
will be (#0, b). 

Let the final configuration be F. Also, S and F agree on the points (-i, -1), (-k, 
-1) and (-k, 0) for all 0 < i < /c — 1. Therefore, hdN{F{—i, —1)) = pq. From the 
definition of Sr, bdw{F{—k,0)) = {#i,c) and the hdw{F{—{k — 1),0)) = {x,y) 
where y depends on the value of x. If 3m < k — 1 where #i = #0, and for all 
n > m, i > 0, then for all m > n we have bdE{F{—n,0)) = {#{i — l),c) and 
for all 0 < n < m, we have bdE{F{—n, 0)) = (#0, b) and v{—n, 0) = (x, c). After 
the relationship is determined and right boundary tile IR is attached at (1, 0), 
the insertion process will terminate. Additionally, we have r = 2 which means 
that only tile will attach at a time; the one which has two neighboring tiles. To 
sum up, this clearly indicates that we will have unique tile configuration which 
takes k 1 steps to finish ■ 

5.2 SQUARE-ROOT ASSEMBLY 

Suppose we wish to find square-root of a number then we write its bits on left 
side in input tiles. If the length of n is odd then we add an extra zero bit in 
the leftmost tile and then start filling the bits. The right bits of the input tile 
system has 0 and 1 in the two leftmost tiles respectively and all other tiles have 
zero in the right part. Fig. p!5] show the input configuration for n = 42.25. In this 
system, there are four operations namely compare (red tiles), shift (orange tiles), 
subtract and shift (blue tiles), and insert bit (green tiles). Compare and Insert 
operations have assembly growth from west to east while the others have it from 
east to west. As shown in Fig.[^ we compare input numbers, p = 101010.01 and 
q= 010000.00, from west to east and since p>q subtract and shift tile system 
attaches. After the growth of this tile system, we insert 1 at 2^^ position from 
left. Now, the process is repeated that is compare, shift or subtract and shift, and 
insert perpetually unless the northern glue of insertion tiles becomes equal to 
$(n-l) and then the termination end tiles self assemble themselves. At any point 
if compare tile give < as output then instead of subtract and shift boundary 
tile, shift boundary tile will attach to CR. The boundary tiles for all these tile 
systems are shown in Fig. 

Note that by rational number we mean that the number should be recurring 
or terminating. However, in the recurring case, like calculator, we take first 8 
bits after decimal. 

Theorem 2 Let be an alphabet defined as 

Y, = {0,1, 00,01,10,11, e, /, t, <, >, =, r, dot) IJ 
[xy I X e{<, >, =}; y e{0,1}} IJ 

{ 000 , 001 , 010 , on, 100 , loi, no, 111 } j 

{$i| — 1 < i < n — 1; i = 2k ^ 1} J 

{Ai|l < i < n — 1; i = 2k ^ 1} J 

{#ib\l < i < n - 1; 6 G {0,1}} J 

{$ix\ — 1 < i < n — 1; X G {0,1, s}} 



Fig. 14: Boundary tiles for compare, shift, insert, and subtract & shift (clock¬ 
wise). Here, the variables used in the glue colors of boundary tiles are : 
{$ix| — l<i<n — 1 ;xG{0, 1, 5 }}, {Xi\l < i < n — 1] i = 2k + 1}, 

{#z&|l < z < n- 1;6 e {0,1}}, {xy\x e {<,>,=};y e {0,1}}, {x\x e {>,=}}• 


and T 2 be the tile set defined over SeeTable 6.1, given that r = 2 and g=l, 


can compute square root of a n bit rational number defined as Sr = (T 2 , S^g^r) 
where S is the seed configuration.. Then, there exists some {xo,yo)eZ 2 such that 
: Sr{xo - k,yo - 1) = So, SR{xo,yo - 1) = Eq, Sr{xo - k,yo) = CL, and for 
all i e{0,1, ...,k — 1}, bdN{xo — i,yo — 1) = pq. Then, the configuration S can 
produce unique final configuration F and will compute square-root of a number. 


Proof. Consider the tile system Sr. Let input tiles be T 2 = Eq = {e, null, null, null}, 
Doo = {00, null, null, null}, D 01 = {01, null, null, null}, 

Dio = {10, null, null, null}. Dll = {11, null, null, null}. So = {$ix, null, null, null}}. 


Let the seed configuration S : Z 2 ^ T 2 be: 

^Eo = Sr{1,-1) 

So = SR{-k,-l) 

< CL = SR{-k,0) (3) 

\/ie{0,1,.., k - 1}, SR{-i, -1) = Dpq 
For all other (x, y)eZ 2 , SR{x,y) = empty 

For the given seed configuration S as shown in Fig.[^ there is only on position 
where the incoming tile can attach. Also, given that r = 2, by induction we 
can say that Vte T 2 , {bds{t),bdw(t)) and (bds{t),bdE{t)) is unique. Therefore, 
it follows that Sr produces unique final configuration E and it takes kFl steps 
to assemble. 

Let the final configuration be F. Also, S and F agree on the points (-i, -1), 
(-k, -1) and (-k, 0) for all 0 < i < /c — 1. Consider the seed configuration as shown 
in Fig. 1^ According to Lemma 1 of [T], we will get unique configuration for the 















































































Fig. 15: Square Root of n = 42.25 is calculated by the self-assembly of tiles. 
Here n is encoded in binary as 101010.OI 2 . We ge 110.1 2 as the output at the 
left side of the tile assembly model. Color of the tile represent same tile-set as 
mentioned. Yellow tiles simply indicates decimal point. Grey tiles indicate values 
after decimal point. 


comparison tile system which will terminate with CR attaching at (1, 0) position. 
Now, if bd^iCR) is = or >, then only SSR will attach to this configuration since 































































































































































































it is the only matching tile. Now, we have new seed configuration which resemble 
to subtract and shift seed configuration. Again, according to Lemma 3 of [1] , this 
will produce unique final configuration in /c +1 steps which terminates when SSL 
tile attaches at (-k, 0). If on the other hand, bd^iCR) is < than SR tile will 
attach. In the case also, according to Lemma 2 of [T], we will get a unique tile 
configuration in /c + 1 steps. This process terminate when SL will attach at (-k, 
0) position. After either of this configuration is formed, IL will attach at (-k, 
-1) position since we have r = 2 which makes it the only possible tile to be 
attached to this growing assembly. According to Theorem 1, this insertion seed 
configuration will produce unique output tile assembly. 

Now, the same process repeats unless we get the desired number of bits after 
decimal place. Therefore, according to induction, the entire configuration will 
be unique. After all the computation is finished, the corner tiles and boundary 
tiles will attach at the top to stop the growth of assembly. Hence, square-root 
computation is terminated ■ 

Theorem 3 The time taken for eomputing square-root is 0[{ki + ^^ 2 )^) where 
ki is number of bits before deeimal and k 2 is number of bits after deeimal. 

Proof. In our previous proofs, we have already shown that all the tile as¬ 
semblies are step configuration which means that only one tile is attached to the 
seed configuration at a time. Also, since r = 2 no tile will get attached to the 
configuration unless it has two surrounding tiles. Also, our tile assembly takes 
kc levels to compute and in each level we have n = ki k 2 computations. This 
makes total computation time to be equal to kc ^ {ki + /C 2 ) where kc = 3n/2. 
Therefore, since n = + /c 2 , we have running time complexity as 0((/ci + /c 2 )^)B 


Basic Tiles 

Tile Types 

Left Frame 

2n 

Right Frame 

11 

Corner 

4 

Top Frame 

4 

Input 

4 

Computational 

12 T 8 T 16 T 477/ 


Table 1: Tile Set to compute square-root of n bit number. The number n is of 
even length and if not then made even. 

6 APPLICATIONS 

The method of square-root can be used to compute rational numbers. Also, a 
slight modification can help us to compute an irrational number tt. 

6.1 COMPUTING RATIONAL NUMBERS 

Suppose we are interested in finding the square-root of 1/3 then, first, we need 
it in decimal form. Therefore, we propose a method to compute rational number 










p/q given two input integers p and q. It is an infinitely growing tile assembly 
which is bounded from west and south, the other ends being open. The method 
to compute rational number requires some new additional computational tiles 
but, fortunately, this method will not require the boundary tiles. Fig. p!7| 
and [T^ shows additional tiles required to make an infinite assembly growth. The 


growth for computing 1/3 is shown in Fig. 19 



Fig. 16: Additional 9 computational tiles required for comparing bits. Here x 
e{<, >, =} and y e{0,1}. 



Fig. 17: Additional 2 computational tiles required for shifting divisor bits by one 
position. 




Fig. 18: Additional 2 computational tiles required for subtracting and shifting 
bits. Here, x e{>, =}. 


Theorem 4 Let ^ be an alphabet defined as 











































E = { 0 , 1 , 00 , 01 , 10 , 11 , T,t, <,>,=} u 

{xy I X e{<, >, =}; y e{T, 00, 01,10,11}} U 

{000,001,010, on, 100, loi, no, 111} 


and T 3 he the tile set defined over See Table 6.T given that r = 2 and 


g=l, ean eompute a rational number from a given numerator and denomina¬ 
tor. Therefore, the tile assembly is defined as defined as Sr = {Ts,S,g,r) where 
S is the seed eonfiguration.. Then, there exists some (xo,^o)e ^2 sueh that : 
Sr{xo - k,yo - 1) = So, SR{xo,yo - 1) = 00, Sr{xo - k,yo) = CL, and for 
all i e{0,1,..., k — 1}, bd]y{xo — i,yo — ^) = pg- Then, the eonfiguration S ean 
produee a unique and unbounded final eonfiguration F and will eompute rational 
number. 



Fig. 19: Here we divide p = 1 by q = 3 to calculate the rational number 0.3333.... 
Here p is encoded in binary as IO 2 and q is encoded as 112 . We get 0.0101 ..2 
as the output at the left side of the tile assembly model. Arrows indicate the 
tile assembly is unbounded while triangle indicates that the entire space will be 
filled with ( 0 , 0 ) tiles with relevant glues. 
























Proof. Consider the tile system Sr. Let input tiles be T = 
{Dqq = {00, null, null, null}, Sr = {null, = ss, null} 

Dqi = {01, null, null, null}, Diq = {10, null, null, null}. 

Dll = {11, null, null, null}. So = {$ix, null, null, null}, 

00 = {T, null, null, null}}. 


Let the seed configuration S : Z 2 ^ T he: 

'00 = Sr{1,-1) 

So = SR{-k,-l) 

< CL = SR{-k,0) (4) 

Viejo, 1 , k - 1}, SR{-i, - 1 ) = Dpq 
For all other (x, y)eZ 2 , SR{x,y) = empty 


Let the final configuration be F. Also, S and F agree on the points (-i, -1), 
(-k, -1) and (-k, 0) for all 0 < i < /c — 1. The seed configuration S for the tile 
growing tile assembly shown in Fig. 13 resembles the original compare, shift, 
and subtract/ shift configurations. There is only on position where the incoming 
tile can attach but, in this case, after the terminating tile has attached we can 
have growth on both of its sides. Left side of the terminating tile is exactly 
similar to division seed configuration hence, according to Lemma 4 of [1], it will 
have unique growth. As far as the growth of right side of terminating tile is 
concerned, there are only zero tiles with all the glue colors having values like 
{000,00,0, < 0, > 0, = 0} depending on their level. 

The same process continues perpetually. Therefore, according to induction, 
the entire configuration will be unique. Therefore, we get an infinite recurring 
number. If the rational number is terminating then we will get trailing zeroes. 


Type of tiles 

Number of Tiles 

Compare 

12 + 9 = 21 

Shift 

8 + 2 = 10 

Subtract & Shift 

16 + 2 = 18 

Left Boundary 

3 


Table 2: Tile Set to compute rational number p/q. 

6.2 COMPUTING IRRATIONAL NUMBER - PI(7r) 

In the final section of the paper, we will demonstrate a way to compute tt using 
an infinite series. The computation method uses same tile systems - Subtraction, 
Shift, Duplicator, Insertion - to compute the value of an infinite series. In order 
to compute the value of tt, we have used Gregory Leibniz series which is defined 
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as: 











To compute this series, we will require values of all the fractious so that we sum 
them up. To hud the value of fractious, therefore, we cau use Theorem 4 which 
cau compute the value of fractiou. Ouce we have the uumber iu decimal form, 
we simply ueed to do additiou or subtractiou to already existiug series sum. 

We ueed the value of output of divisiou iu simpler form so that we cau use 
it for further computatiou. To do this, we have used copy aud shift tile systems. 
Refer Fig. The value of quotieut ou the left bouudary of divisiou tile set is 
extracted by shiftiug them oue bit per level - yellow aud orauge tiles - as showu 
iu Fig. As a result, wheu we fiuish the divisiou process, we have the output 
value aloug with existiug sum to which we either waut to add or subtract this 
fractiou. Now, depeudiug ou whether the term is odd or eveu we ueed to add 
or subtract it to the existiug sum which is showu by greeu tiles iu Fig. This 
arithmetic occurs meauwhile we have iuput tiles - grey colored - for uext divisiou 
beiug attached ou the other side. Ouce this additiou or subtractiou is fiuished, 
we replace the curreut sum with this uew sum iu the grey tile as showu iu Fig.|^ 
This process keeps ruu iu perpetuity. Depeudiug ou the uumber of precisiou bits 
we cau coutrol the growth of our assembly. 

Theorem 5 Let be an alphabet defined as 

E = {0,1,00,01,10,11, T, t, <, >, =, +, -} u 

{xy I X e{<, >, =}; y e{T, 00 , 01 , 10 , 11 }} U 

{000,001,010, on, 100, loi, no, 111} u 

<i <n-l; i = 2k + l} [J 

{#ix\ — 1 < i < n — 1; X G {0,1, s, ss}} 

and T 4 be the tile set defined over The tile set ineludes eompare, insert, shift, 
and subtraet & shift tile systems. Given that r = 2 and g=l, ean eompute an 
irrational number n. Therefore, the tile assembly is defined as defined as Sr = 
(T4, S,g,r) where S is the seed eonfiguration.. Then, there exists some {xo,yo)eZ 2 
sueh that : SR{xo,yo) = So, Sr{xo, + 1 ) = CL, and for all i e{ 0 ,1,..., k — 1}, 
bd^ixo — i,yo — = PQ- Then, the eonfiguration S ean produee a unique and 

unbounded final eonfiguration F and will eompute irrational number - tt. 

Proof. Siuce we have used the same tile systems - compare, shift, subtract 
aud shift, aud duplicator - its proof will be similar to previous theorems. There¬ 
fore, we will uot go iuto details of the proof ■ 

There are a few observatious about GregoryLeibniz series [20] . 

— The series couverges at a very slow rate [ 20 ]; Leibuiz’s formula couverges 
extremely slowly: it exhibits subliuear couvergeuce. Calculatiug tt to 10 cor¬ 
rect decimal places usiug direct summatiou of the series requires about five 
billiou terms. 

— The series is easy to compute. The uumerator of each term is 4 while the 
deuomiuator iucreases by 2 at every-step. 

— The series does uot require multiplicatiou which meaus we do uot require 
auy uew tile system. 


00 3r[; _I0_ 3rC 10 w: dot w: 10 3 t-C 00 3K 00 3rC 10 3rC 00 3rt: 00 
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Fig. 20: Computing value of an irrational number tt - 3.141592.... Here, grey tiles 
are input tiles, green are addition/subtraction tiles, yellow are copy/duplicate 
tiles, and orange tiles are shift tiles. Blue tiles indicate the reminder of the 
corresponding division process. 




























































































































































































Note that the method to compute tt may not be most optimized version. As a 
result, it is an open problem to compute value of tt in an optimized fashion. 

7 SUPPLEMENTARY FILE 

We have also written a program to generate .tile file of square-root which can 
be downloaded from our website http://www.guptalab.org, 

8 CONCLUSION 

Our method to compute square-root, and rational numbers is simply a modifi¬ 
cation of the Xuncai et. al. method to find subtraction and division, using DNA 
tiles. Our system uses 0(n) computational tiles for square-root of an n bit bi¬ 
nary number and 0(l) tiles for computing rational number. However, it would 
be an interesting to compute tt using finite number of tiles. 
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